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ABSTRACT 

Collisionless shocks are commonly argued to be the sites of cosmic ray (CR) acceleration. We study the influ- 
ence of CRs on weakly magnetized relativistic collisionless shocks and apply our results to external shocks in 
gamma-ray burst (GRB) afterglows. The common view is that the transverse Weibel instability (TWI) generates 
a small-scale magnetic field that facilitates collisional coupling and thermalization in the shock transition. The 
TWI field is expected to decay rapidly, over a finite number of proton plasma skin depths from the transition. 
However, the synchrotron emission in GRB afterglows suggests that a strong and persistent magnetic field is 
present in the plasma that crosses the shock; the origin of this field is a key open question. Here we suggest that 
the common picture involving TWI demands revision. Namely, the CRs drive turbulence in the shock upstream 
on scales much larger than the skin depth. This turbulence generates a large-scale magnetic field that quenches 
TWI and produces a magnetized shock. The new field efficiently confines CRs and enhances the acceleration 
efficiency. The CRs modify the shocks in GRB afterglows at least while they remain relativistic. The origin of 
the magnetic field that gives rise to the synchrotron emission is plausibly in the CR-driven turbulence. We do 
not expect ultrahigh energy cosmic ray production in external GRB shocks. 

Subject headings: acceleration of particles — cosmic rays — gamma-rays: bursts — plasmas — shock waves 



1. INTRODUCTION 

Collisionless shocks are observed on all astrophysical 
scales. The diffusive shock acceleration (DSA) mechanism 
is believed to accelerate cosmic rays (CRs) in these shocks 
dBelji ri978: Bla ndford & Ostrikerlfl978l: IBlandford & Eichleri 
fl987h . The CRs can carry a substantial fraction of the en- 
ergy of the shock and thus t he CR pressure can influence 
the structure of the shock dEichlerl 1 1 979t IBlandford! 119801: 
lDrurv& Vol5[l98Tt lEllison i& Eichleri 1 19841) . This picture 
has recently received sup port from X-ray ob servations o f su- 
pernova (SN ) remnants dWarren et al.ll2005l) . Recently. iBelll 
(12004 120051) has shown that CRs in SNe can drive turbu- 
lence and amplify magn etic fields in the shock upstream, and 
iDar & de Rujulal (120051) speculate that relativistic jets could 
do the same. 

The current lore is that weakly magnetized relativistic col- 
lisionless sh ocks are mediated by the t ransverse Weibel insta- 
bility (TWI: IWeibell959l: iFriedll 1 9591) . TWI produces a mag- 
netic field near equipartition that provides collisional coupling 
in the shock transition layer dGruzinov & Waxmanl 119991 : 
iMedvedev & Loebl H999). TWI is universally observed in 
two- and three-dimensional particle-in-c ell (PIC) simulations 
of unmagnetize d colliding shells ( e .g.. iLee & Lampel 1 1973 1: 
Gruzinov 1 12001b ISilya etaD l2003h iFrederiksen et all 120041 
Jaroschek et al. 2004: iNishikawa et alJ l2005: Med vedev et alJ 
20051: iKatol 120051) . The resulting field is small-scale, of the 



order of the proton plasma skin depth A s . 

The role of TWI in particle acceleration and in the down- 
stream thermodynamics is controversial. The field is expected 
to decay rapidly , within a few A s from the sh ock transition 
(lGruzinovll20"0U iMilosavljevic & Nakaril2006l) . This decay 
is evi dent in three dimensional PIC simulations ( Spitkovsky 
120051) . Although there are claims that th e decay saturates 
at distances of — 100A S from the shock dSilva et al.1 120031: 

' Theoretical Astrophysics, Mail Code 130-33, California Institute of 
Technology, 1200 East California Boulevard, Pasadena, CA 91 125. 
2 Hubble Fellow. 



IMedvedev et al.ll2005l) . survival of the field over larger dis- 
tances has not been demonstrated. 

According to the pop ular model of gamma-ray bursts 
(GRBs; e.g., Piran 2005, and references therein), the after- 
glow originates in a relativistic blast wave that propagates into 
an ambient e~p plasma. The afterglow emission is ascribed 
to synchrotron radiation from nonthermal electrons that gy- 
rate in the shock-generated magne tic field. Detailed studies 
of GRB spectra and light curves (Panaitescu & Kumari 120021 : 
I Yost et al]2 003) have shown that the magnetic energy density 
in the emitting region (the downstream shocked plasma) is a 
fraction es ~ 10~ 2 to 10~ 3 of the internal energy density. 1 This 
field must persi st over at least a few percent of the width of 
the blast wave dRossi & Reesll2003l) . Its origin has remained 
a key open question. 

Compressional amplification of the weak pre-existing mag- 
netic field of the interstellar medium (ISM) yields eg ~ 10~ 9 
dGruzinov! 1200 ll) and does not explain the field in the emit- 
ting region. If TWI does develop in the transition layer, it 
generates a strong field in the vicinity of the shock transition. 
However, as discussed above, there is no evidence that it can 
persis t over the required ~ 10 9 A S from the shock (e.g.. lPiranl 
2005 and refer ences therein). 

X-ray observations of GRB afterglows are modeled as 
an optically thin synchrotron spectrum requiring nonthermal 
electrons with Lorentz factors as high as ~ 10 6 . The obser- 
vations indicate that electrons, and thus protons as well, are 
efficiently accelerated in the shock to produce a hard power- 
law spectrum. DSA can achieve such Lorentz factors if the 
circum-burst medium is magnetized at ~ 1 /iG level, as ex- 
pected in the ISM. The weak magnetic field, however, does 
not directly affect the structure of the shock (e.g., TWI can 
still develop in the shock transition). 

It seems that a self-consistent picture of relativistic astro- 
physical shocks includes high energy particles (CRs) and at 

1 Under some circum stances tg as low as ~ 10~ 5 can fit the data 
(Eichler & Waxman 2005). 
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least a weakly magnetized upstream. Here we explore the in- 
teraction between these components. In § [2] we derive the 
conditions under which the CRs drive turbulence in the up- 
stream; this turbulence generates a large-scale magnetic field 
that increases the acceleration efficiency. 2 In § [3] we argue 
that this mechanism is a candidate solution of the origin of 
the field inferred from the afterglow emission in external GRB 
shocks. In § |U we discuss implications for TWI and for the 
production of ultrahigh energy cosmic rays (UHECRs). 

2. COSMIC RAY - MAGNETIC FIELD INTERACTION 
2.1. Cosmic Ray Trajectories and the Return Current 

Consider a relativistic shock wave with Lorentz factor Y 
that accelerates particles to Lorentz factors 7 ^ T. The ac- 
celerated particles are roughly isotropic in the downstream 
frame; a fraction of them are moving ahead of the shock. In 
the upstream frame, the velocities of the accelerated particles 
are directed within an angle 6 ~ V2/Y of the shock normal. 3 
If the upstream contains a magnetic field B, the field deflects 
charged particles moving ahead of the shock. When the veloc- 
ity of a particle is deflected to an angle ~ 1 /Y from the shock 
normal, the shock catches up and re-absorbs the particle. 

The orbit of a particle in the upstream field can be quasi- 
circular or diffusive, depending on how the radial distance 
X the particle traverses ahead of the shock compares to the 
relevant magnetic correlation length A. The typical angle by 
which the magnetic field deflects the particle is AO ~ X / r g 
when the orbit is quasi-circular and Ad ~ (AX) 1 / 2 /r g when 
the orbit is diffusive, where 



seed field 



jmc 
eB(\) 



(1) 



is the gyroradius of the particle expressed in terms of the mag- 
netic field B(X) on scales A. Reabsorption in the shock occurs 
when A8 ~ 1 /Y, and thus the distance traversed by the parti- 
cle in the quasi-circular and diffusive regime is 



X 



(r g /D 2 /A, (X»A) 



(2) 



In spherical blast waves, a particle can be accelerated by re- 
peated shock crossing if the time separating its emission and 
re-absorption in the shock ~ X/c is smaller than the expan- 
sion time of the shock ~ R/c, where R is the shock radius. 4 
Therefore X < R is required for acceleration. The maximum 
distance that the particle reaches from the moving shock tran- 
sition is substantially smaller, 



A 



X 

f2- 



(3) 



because the particle and the shock are both moving at about 
the speed of light. Another requirement for acceleration is 
that the cooling time be longer than the acceleration time (see 
§0. 

In general, protons and electrons are both emitted into the 
shock upstream. However, the synchrotron and self-Compton 
cooling inhibit the acceleration of electrons, and so the pro- 
tons, which we refer to as the CR precursor of the shock, can 



2 IGruzinov (2001) and Waxman 12004) speculate that CRs could amplify 
large-scale magnetic fields in external GRB shocks. 

3 Unless otherwise noted, quantities are evaluated in the upstream frame. 

4 It is also implicitly assumed that some fraction of accelerated particles 
can diffuse back to the shock after they find themselves in the downstream. 
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FIG. 1 . — Schematic representation of the reaction of the upstream plasma 
to the CR streaming ahead of the shock. The figure is in the frame of the 
shock so that the upstream plasma moves vertically downward. An initially 
weak magnetic field loop is stretched by the Ampere force ( J re i X B)/c and 
plasma is accelerated sideways as it flows toward the shock front. 



reach higher energies and farther from the shock into the up- 
stream. The upstream fluid sees a net positive CR current 
|Jcr| ~ <?v s «cr in the region where it overlaps with the pre- 
cursor. Here v s and «cr are the CR velocity and density. In 
the upstream region, which is not accelerated to relativistic 
velocities by the CR pressure, the CR velocity is about the 
speed of light, v s ~ c; thus, we have 



Jcr ~ ec« C R. 



(4) 



To maintain overall charge neutrality, the upstream re- 
sponds to the CR current by supporting an opposite "return 
current" J ret that short-circuits the electric field associated 
with the charge separation in the CR precursor. 5 If «cr is 
smaller than the upstream density, the return current approxi- 
mately balances the CR current, 



' Jcr- 



(5) 



If the upstream contains a weak magnetic field, the return cur- 
rent couples with this field; next we explore the consequences 
of this coupling. 

2.2. The Driving of Turbulence and Field Amplification 

Here we argue that the return current drives turbulence in 
the upstream and suggest that this turbulence leads to mag- 
netic field growth on scales much larger than the proton 
plasma skin depth. One type of interaction between the CR 
precursor and the magnetic field of the upstream, explored by 
iBelll d2004 12005b in the context of Newtonian shocks, is the 
Ampere force 

F=^L^ (6) 

c 

which accelerates the upstream perpendicular to the shock ve- 
locity. We assume that the initial weak magnetic field has 
power on all scales, as expected in interstellar turbulence, and 
that the power on relevant scales is roughly scale-independent. 

Consider a loop of radius A of a weak initial magnetic field 
Bo parallel to the shock transition that is directed clockwise as 
seen from the shock, as shown in FigureQ] and assume that the 
upstream is initially stationary. The Amper e force accelerates 
the fluid away from the center of the loop. IBelll {2004, 2005) 
carried out magnetohydrodynamic (MHD) simulations of this 
process; we base our estimates on his results. 

5 The return current is established on the plasma time < 10 ms which is 
instantaneous compared to the light crossing time of the CR precursor. 
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The dynamics of the upstream fluid can be approximated 
using MHD, and thus flux freezing during the expansion of 
the loop implies that 



Be 
pr 



= const, 



(7) 



where Bg and p are the azimuthal magnetic field and the 
density at radius r, respectively. If the expansion is non- 
relativistic and if the thermal and magnetic pressure are ig- 
nored so that only the Ampere force influences the motion, 
the radius of the loop accelerates according to 



d 2 r 
dt 2 



JcrBo 



(8) 



pc Xpoc 

Therefore, the loop expands exponentially r(t ) <~ Xe a ' at a rate 



/ JcrBu 
V Xpoc 



1/2 



(9) 



implying a velocity of expansion of dr/dt ~ ra. 

The upstream accelerates until the pressure overtakes the 
Ampere force (we assume that the expansion remains sublu- 
minal at this point). In a realistic environment where the mag- 
netic field fluctuates on all scales, we expect the expansion on 
each scale to saturate when neighboring magnetic "rings" on 
similar scales collide, namely when r ~ 2A. This happens 
when 

crA 

at 1. (10) 

c 

Compressional amplification of the magnetic field at this point 
is of the order of unity. However it is commonly argued 
(Kulsrud 2005, and references therein) that a turbulent dy- 
namo produces magnetic energy near equipartition with the 
turbulent energy. Therefore, the upstream fluid that is exposed 
to a current Jqr over a time t will develop turbulent motion on 
all scales 

A <A„« (11) 
poc 

Each eddy amplifies the magnetic energy on its own scale 
on an e-folding time l/cr(A). The bulk of the kinetic en- 
ergy initially stored in expanding shells ~ j(crX mRX ) 2 turns 
into turbulent motion and the magnetic field on scales ~ A max . 
The field in equipartition with the energy in turbulent eddies 
equals 

/4,AW /2 
V cpo J 

This is the minimum field strength generated in the shock up- 
stream on scales A < A max . The interaction between the CRs 
and the generated field could continue to accelerate the fluid 
after the shells have collided once. This could result in ad- 
ditional turbulent driving and field growth. Then a saturation 
occurs when the Ampere force balances the magnetic tension 



|F|~ — 

which implies a limit on the final field of 

47rA7 re t 4vrA7cR 

B 2 ' 



(13) 



(14) 



c c 

iBelll (|2004), who considered Newtonian shocks, argues that 
such a field is generated as a consequence of CR streaming. 



Therefore we can limit the generated magnetic field to lie be- 
tween B\ and B 2 . 

For some A, the above values of B formally imply superlu- 
minal motion. We do not analyze the evolution of the mag- 
netic field after its energy approaches equipartition with the 
rest energy of the upstream. 6 Our treatment is applicable at 
distances from the shock where the magnetic field does not 
reach equipartition, B < (47rp; sm ) 1//2 c. 

2.3. The Cosmic Ray Spectrum and the Quasi-Steady State 

The generated magnetic field can modify the spatial pro- 
file of 7cr by confining CRs closer to the shock than the pre- 
existing field. The maximum coherence length of the new 
field cannot exceed the maximum distance CRs reach from 
the shock ~ R/T 2 . Therefore the propagation of the most en- 
ergetic CRs traversing a distance X ~ R between emission and 
re-absorption is always diffusive in the new field. If the CR 
motion is dominated by the new field, 7 the shock settles in a 
quasi-steady state in which to every distance A < R/T 2 from 
the shock corresponds a 7 such that CRs with Lorentz fac- 
tor 7 are confined within a distance A from the shock by the 
magnetic field generated by the streaming of these same CRs. 

The density of the CRs at a distance A is 

" Cr(7) ~4^A7t^' (15) 

where dN/dj is the CR energy spectrum. The CR spectrum 
produced by the DS A is typically a power law 

dN 



c/7 



oc 7 



(16) 



< -7 < 



7 — 7max- 



We take the minimum Lorentz 



in a range y n , 

factor of CRs to equal ( Achte rberg et al.|[200lL and references 
therein) 

7min~r 2 . (17) 

Numeric al simulations of acce l eration in ultrarelativis- 
tic shocks dAchterbergetal.1 120011: lEllison & Double! 120021: 
Lemo ine & Pellefieri 120031) agree with the spectral index 
p_=_J2 Rj 2.22 derived analytically by iKeshet & Waxmanl 
(2005) for isotropic diffusion. While the true spec- 
tral index will depend on the detaile d interaction be- 
tween CRs and magnetic turbulence (e.g ., Elliso n~& Double! 
2004tlLemoine & Revenull2006l: iNiemiec & Ostrowskill2006l: 



Lemoine, Pellet ierT Revenul 120061) and may differ from this 



derived value, we adopt 



20 



(18) 



in what follows. If Ecr is the total energy in CRs in the shock 
wave, then the CR spectrum is normalized such that 



Hi 



"fnipC c/7 

dj 



CR- 



(19) 



We relate the L orentz factor of the shock to the total energy 
in the blast wave (Blandf ord & McKedl 19761) 



r : 



: 2 Ri 



L/2 



(20) 



6 Magnetic field exceeding equipartition with the rest energy of the up- 
stream would affect the hydrodynamic profile of the shock and accelerate the 
upstream in the direction of shock propagation, thereby reducing v s and 7cr- 

7 The preexisting magnetic field could have power on the largest scales ~ R 
and thus it could dominate the CR motion even if the new field is stronger. 
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where no « po/m p . We assume further that a fraction 



Ecr 
E tol 



(21) 



of the blast wave is stored in the CRs. 

We proceed to outline the procedure by which the shock- 
accelerated CR Lorentz factor, the shock-generated magnetic 
field coherence length, and the energy density in the shock 
generated magnetic field are calculated self-consistently as a 
function of distance from the shock A. In §[3] we apply this 
procedure to external GRB shocks. 

2.4. The Self-Consistent Solution 

At any distance A from the shock, the CR current is dom- 
inated by the least energetic CRs that reach this distance. 
CRs with smaller 7 are confined to shorter distances from 
the shock. The CR current 7cr(A) is larger at smaller A. 
Note the time over which an upstream element is exposed to 
7cr(A) is proportional to A <C R, i.e., since the cosmic rays 
are confined so close to the shock front in an ultrarelativis- 
tic shock, the time available for generation of the magnetic 
field is short. The field is generated on scales A for which 
a(X,A)A/c > l.. 8 In calculating a, the distance CRs reach 
from the shock and the CR number and current densities are 
all calculated using the generated field B, rather than the pre- 
existing field Bo. A self-consistent solution for the CR and 
magnetic field profile ahead of the shock shows that CRs with 
smaller 7 generate a stronger magnetic field (because of their 
larger Jcr) on smaller scales A (because of the shorter expo- 
sure time) and are confined by their self-generated field. The 
maximum CR Lorentz factor j m - dx and the maximum field cor- 
relation length A max are obtained for A ~ R/T 2 (see § 12. j} . 

To solve for 7(A), A(A), and B(A), we proceed as follows. 
First, we eliminate r g and X from equations (Q~|i, (f2j the diffu- 
sive case), and (0 to obtain 



AA 



2 2 4 

7 m c 



(22) 



T 4 e 2 B 2 

Next, we eliminate «cr, 7min, and Ecr from equations and 
dBJ-EB to obtain 

r 4 / 9 e C R£,ote 



JCR ■ 



l87rj n / 9 m v cR 2 A' 
Next, we eliminate a from equations (O and ( TTOb to obtain 



(23) 



/JcrB \ 1/2 A i 
V Apoc / c 



(24) 



Finally, we set the strength of the saturated magnetic field to 
B2 where the Ampere force balances magnetic tension (see 
§ 12.21 the results for saturation at B\ are qualitatively the 

same) 

B (25) 
c 

We substitute T from equation (|20T > in equations (l22l and 
( f23l . Then we substitute Jqr from equation d23"b in equations 
d24l > and 1251 . We finally solve for 7, A, and B as a function 
of A and the parameters of the blast wave. It is convenient to 
express the magnetic field strength in terms of 



(B 



4irpc 2 



(26) 



and to express distance from the shock in terms of the dimen- 
sionless parameter 



. Ar 2 

A= < 1. 

R ~ 



(27) 



The maximum 7 and A correspond to CRs that reach farthest 
from the shock, i.e., to A = 1. 

3. GAMMA-RAY BURST AFTERGLOWS 

External shocks in GRBs are the best astrophysical candi- 
dates of weakly magnetized relativistic collisionless shocks. 
They have Lorentz facto rs r g", 10 when the blast wave is 
at a radius R ~ 10 17 cm (lPiranll2005l and references therein). 
The blast wave decelerates and becomes Newtonian at R ~ 
10 1S cm. The external shocks are initially beamed in the form 
of a jet with a typical opening angle 6} ~ 0.1 rad. The angular 
size of the causally connected region in the shock is 1 /V. At 
early times when T > 1/6}, the blast wave has isotopic equiv- 
alent energy E tot = 10 52 - 10 54 ergs and evolves as a spherical 
fragment with a Lorentz factor given by equation ( 1201 ). The 
evolution at late times when T < 1/6} is poorly understood: 
the opening angle is expected to increase; T decays faster 
(perhaps exponentially) with R, and the isotropic equivalent 
energy approaches ~ 10 51 ergs. 

Here we explore the interaction between the CRs and the 
pre-existing magnetic field in the circum-burst medium while 
the evolution of the blast-wave is quasi-spherical (V > 1/6}). 
We expect that our results are qualitatively applicable also 
when r < 1/6} by taking E tot - 10 51 ergs and R ~ 10 18 cm. 
We derive the maximum Lorentz factor 7 of CRs that reach a 
distance A from the shock, the maximum magnetic field cor- 
relation length A at this distance, and the fraction eg of the 
energy density in the magnetic fields on sc ales A. 

Following the solution outlined in § !2.4l we obtain 



7' 
A' 



• 5 x 10 7 A 0,25 e°J 2 B°f E° 53 75 R^i n^ 31 



10' 



,0.25 



,0.55 



-0.7 
J 53 



? 3.04 
M7.5 



eg ~ 10 A" 



-1.21 ,0.99 D 0.19 ,71.21 d-3.84 - 



1.19 



M7.5 



(28) 



' In reality, field saturation will require multiple e-foldings, a A/c ~ few. 



where A = 1 corresponds to the highest energy CRs that are 
accelerated in the shock. Here, ccr = 0.1 e_i is the frac- 
tion of the energy E tot carried by CRs ahead of the shock, 
E to t = 10 53 £ , 53 ergs, Bo = B-6 mG, «; sm = «o cm" 3 , and R = 
10 17 5 7?i7 5 cm. We have used p = y but the dependence on 
p is weak. Evidently, in the shocks considered here, the field 
generated in the precursor has the capacity to confine CRs 
with higher energies than a preexisting microgauss field. 9 

7 —1 /2 

Note that A max > A s ~ 2 x 10 n ' cm. The fractional 
magnetic energy density e# should be approximately pre- 
served as the fluid passes the shock transition. The depen- 
dence of eg on the radius and the field correlation length is 
€b (x R lAl X~ 113 (it is independent of E tot if ecR is constant). 
Note that eg is larger on smaller scales since it is generated 
at smaller A by the larger 7cr that is dominated by CRs with 
smaller 7. 

Since the upstream is turbulent on scales < A max , and since 
this turbulence is accompanied by density inhomogeneities, 
the shock transition and the fluid passing the transition are 
also expected to be turbulent on the same scales. The evo- 
lution of the magnetic field after it passes the transition is 

9 The maximum Lorentz factor to which CRs are accelerated in the pre- 
existing field of the ISM is j msx .(B ) ~ 4.6 X 10 6 fi]' 2 B_ 6 n Q [/2 R^. 
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described by the laws governing MHD turbulence, and we 
do not attempt to describe it here. The value of es in the 
downstream region will be that reached in the shock upstream, 
modified by any evolution that the turbulence experiences as 
it travels downstream past the shock transition. 

The main premise of our treatment is that the precursor is 
dominated by protons and thus carries a net positive current. 
To compete with the protons in the distance they reach ahead 
of the shock, the electrons would have to have Lorentz factors 
7e > 7min'«p/'«e ~ T 2 m p /m e . Electrons with 7 e ~ 10 8 - 10 1() in 
the rest frame of the upstream are cooled by the synchrotron 
mechanism on timescales much shorter than the age of the 
shock, and thus their acceleration to even higher energies is 
suppressed. Therefore, indeed, the shock precursor at proton 
Lorentz factors 7 > 10 5 is completely dominated by protons 
and carries a net positive current. 

4. DISCUSSION 

The picture presented here casts doubt on the role of TWI 
in collisionless shocks. It is commonly argued that TWI gen- 
erates the magnetic field that facilitates collisional coupling 
in the shock transition. However, if the upstream is weakly 
magnetized, the shock accelerates CRs, and these CRs drive 
upstream turbulence and magnetic field generation. When the 
upstream fluid reaches the transition — where TWI is assumed 
to take place — the large-sc ale field generated in the pre cursor 
will plausibly quench TWI. lHededal & Nishikawal(l2005l) find 
that TWI is quenched for ratios of the electron plasma to cy- 
clotron frequency ui pe /ui ce < 5. For nonrelativistic upstream 
electrons this implies quenching for eg > 2 x 10~ 5 , which is 
expected from our analysis. 10 

Recently, a question has been raised whether external 
GRB shocks are sites of UHECR acceleration (7 ~ 10 10 - 
10 11 , iDermerl 120021: IWaxmanl ^OOlh. T he idea that UHE- 
CRs are produced in GRBs dWaxmanl 119951: IVietril 11995b 
Milgrom & UsovjGHH) is motivated by the observation that 
the cosmic energy densities in UHECRs and GRB fireballs 
are comparable. The presumed absence of strong upstream 
magnetic fields in external GRB shocks has focused attention 
on UH ECR production in shocks propagating into relativistic 
eiecta dWaxmanll2004t iGialis & Pelletierll2005l) . Our analysis 
of external shocks indicates that the maximum Lorentz fac- 
tor of CRs produced there is well below the UHECR range. 

10 Because of the finite transverse velocities of the CRs c/Y in the 
frame of the ISM), the CRs themselves will not excite TWI in the upstream. 



Even under the most optimistic assumptions (e^ ~ 1 on scales 
A ~ R/T 2 ), the resulting Lorentz factors are 

7 < 3 x lO 10 ^ 3 ^ 6 , (29) 

where E lot = lO^-Etot.si ergs is the total, beaming corrected en- 
ergy of the blast wave. Therefore we do not expect UHECR 
acceleration in external GRB shocks. 

5. CONCLUSIONS 

CRs accelerated in relativistic collisionless shocks excite 
large-scales turbulence and magnetic field generation in the 
shock upstream. The field generated in the shock precursor 
has power on scales much larger than the proton plasma skin 
depth. The propagation of CRs in the generated field is diffu- 
sive. In external GRB shocks, CRs are accelerated to higher 
energies in the generated field than in the pre-existing field 
of the ISM. The generated field reaches equipartition with the 
energy density in the fluid. The shock transition is turbulent 
with a hydrodynamic profile dominated by CR pressure. The 
commonly invoked TWI is probably quenched in relativistic 
collisionless shocks by the magnetic field and the turbulence 
generated in the shock precursor. External GRB shocks do not 
accelerate UHECRs. PIC simulations of collisionless shocks 
must include a weak upstream magnetic field and simulate a 
spatial domain as large as the CR gyroradius in the upstream 
to observe the acceleration of particles beyond equipartition. 
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